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"o 

^D ' Abstract 

o 

^N) ' The central two-pion exchange A^A'' potential at large distances is studied in the framework of rela- 

tivistic chiral symmetry and related directly to the nucleon scalar form factor, which describes the mass 
flj ' density of its pion cloud. This relationship is well supported by phenomenology and allows the depen- 

^O I dence of the asymptotic potential on the nucleon mass to be assessed. Results in the heavy baryon limit 

are about 25% larger than those corresponding to the empirical nucleon mass in the region of physical 
interest. This indicates that it is very important to keep this mass finite in precise descriptions of the 
A'' A'' system and supports the efficacy of the relativistic chiral framework. One also estimates the contri- 
bution of subleading effects and presents a simple discussions of the role of the quark condensate in this 
problem. 

o' 
o 
a^ 
o 

§ ! 1 Introduction 

QCD is nowadays the main theoretical framework for understanding hadronic processes, but its non-Abelian 
O I character makes low energy calculations unfeasible. The usual strategy for overcoming this difficulty consists 

■ ' in working with effective theories, constructed in such a way as to include, as much as possible, the main 



a 



X 



features of QCD. As in Nuclear Physics most processes involve only the quarks u and d, one requires these 
theories to be Poincare invariant and to possess approximate SU{2) x SU{2) chiral symmetry. The breaking 
of the latter, due to the quark masses, is parametrized by the small pion mass {fx) at the effective level. 

In the last twenty years, effective theories incorporating chiral symmetry have been successfully applied 
to hadronic interactions, including or not electro-weak probes. As far as the hadronic sector is concerned 
it is useful to distinguish two classes of processes, involving only mesons or both mesons and baryons. In 
the case of purely pionic systems, effective Lagrangians are treated relativistically and yield well defined 
power counting procedures, in both pion mass and momenta W, H]. This gives meaning to the idea of chiral 
perturbation theory. When nucleons are present, things become more complicated [^ and calculations were 
performed in the simplified framework of the so called heavy baryon chiral perturbation theory (HBChPT), 
in which nucleons are treated non-relativistically H, 0] . Only recently a well defined power counting scheme 
was proposed for the relativistic approach H. 

In the case of two-nucleon systems, interactions involve processes with a marked spatial hierarchy, whose 
dynamical complexity increases rapidly when one moves inward. In the best semi-phenomenological poten- 
tials existing at present, that can reproduce low-energy observables accurately, the interaction is determined 
by the undisputed one pion exchange, supplemented by a theoretical two-pion exchange potential (TPEP) 
and parametrized at short distances [0, g]. 
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The TPEP is the locus of chiral symmetry in this problem, since it is closely related to the pion-nucleon 
(ttN) amplitude . After the works of Weinberg in the early nineties restating the role of chiral symmetry 
in nuclear interactions M, the interest in applications was renewed and several authors have tackled the 
problem of constructing the TPEP. Initially, non-linear Lagrangians containing only pions and nucleons 
were employed g. However, these Lagrangians do not suffice to describe experimental nN data jl^l and 
the corresponding potential misses even the scalar-isoscalar medium range NN attraction. In a later stage, 
other degrees of freedom were introduced, coherence with ttN information was restored and descriptions 
compatible with asymptotic NN scattering data were produced 0, 11^, O, 111, 11^ . 

It is possible to distinguish two partially overlapping theoretical frameworks in these chiral calculations. 
One of them, adopted in rcfs.[|ll| and [|l3|, is based on HBChPT. In this theory, a rule exists for count- 
ing powers of the typical three-momenta exchanged between nucleons, assumed to be of the order of the 
pion mass. This allows the construction of non-relativistic effective Lagrangians, that include unknown 
couterterms and are used to derive amplitudes in which contributions from loops and these counterterms 
add up coherently. A rather puzzling aspect of these counting schemes, in the case of NN interactions, is 
that they predict a leading contribution in the l/m expansion with the following spin-isospin structure |13|| : 
Vieading — Vs cr^^'' ■cr^^'' + Vr Si2 + Wc t'^'' 't''^-', whcrc T, <T and 5*12 are the usual isospin, spin and tensor 
operators. This means that, in HBChPT, the scalar-isoscalar central potential, which is by far the most 
important phenomenological component of the TPEP, corresponds to just a subleading effect. 

The other approach to the TPEP is typically relativistic |1^ and emphasis is given to the tail of the 
interaction, which is determined by the well known ttN amplitude. No attention is paid to counterterms in the 
effective Lagrangian, that correspond to zero-range interactions and represent processes such as multimeson 
or quark exchanges. Tail and conterterm contributions do not overlap in space and hence can be treated 
separately in NN potentials. 

For systems containing just a single nucleon, relativistic and HBChPT calculations were compared and 
found to be consistent, provided the dimensional regularization scale is set equal to the nucleon mass Km. 
In this work we repeat this kind of comparison, for the leading contribution to the two-pion exchange 
NN potential, and show that relativistic and HBChPT predictions fail to agree by 30%. The reason for 
this important discrepancy may be traced back to the fact that the intermediate ttN amplitude cannot be 
represented by a series around the point that determines the potential at large distances, as shown by Becher 
and Leutwyler p|. 



2 TPEP 

The construction of the TPEP was discussed in detail in ref. ||l^ and here we just sketch the main steps. The 
potential is based on the on-shell NN scattering amplitude containing two intermediate pions, from which 
one subtracts the iterated OPEP, in order to avoid double counting. The isoscalar component, represented 
by T^ , is given by 

2! 7 (27r)4 [fc' - m'] ['fc" - m'] ^^ ^ ' ^^ 

where T^ is the isospin symmetric TriV scattering amplitude for pions with momenta k and fc', Q = (k' + k)/2, 
and the factor 1/2! accounts for the symmetry under the exchange of the intermediate pions |]12[ |8). In 
the sequence we represent initial and final nucleon momenta by p and p' , their mass by m and also use the 
variables q — {k' — k), t — q^ , V — [p' +p)/2m and v — V-Q. 



In general, the amplitude T^ can depend on four independent variables, v, t, k^ and fc'^. As the exchanged 
pions are off shell, one should, in principle, keep fc^ and fc'^ unconstrained everywhere. However, when these 
factors appear in the numerator of amplitudes, we may write fc^ = /i^ + (fc^ — /i^) or fc' ^ = /i^ + (fc' ^ — /i^) and 
use the terms within parentheses to cancel pion propagators in eq.(nh. This kind of cancellation is associated 
with short range terms, which do not contribute to the asymptotic potential. Therefore, whenever possible, 
we neglect them and replace k"^ and k'^ by jj? in the numerators of our expressions. Concerning the variables 
V and t, the conditions of integration in eq.(|l|) are such that the main contributions come from the unphysical 
region i^ « 0, i > 4/i^ [R. In particular, the vicinity of the point t = 4/i^ determines the potential at very 
large distances. 

In order to understand the structure of the TPEP, it is convenient to isolate in r+ a term T^, involving 
only chiral pion-nucleon interactions at tree level and written as 

1^ I m m 1 

(2) 
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r+ = ^ u <! 1 

m 



{p + k)^ — iri^ (p — k/y 



where g is the ttN coupling constant. The remainder is represented by T^ and encompasses other degrees 
of freedom, either in the form of low energy constants or by means of specific models, that may include 
the A and other states. One writes symbolically [T+J — [T^] + [T^] for each nucleon and the potential 
is then proportional to [T+](i)[T+](2) = [T+](i)[T+](2) + {[T+](i)[T+](2) + [r+](i)[r+](2)} + [r+](i)[r+](2). 
This product involves box, crossed box, triangle and bubble diagrams, that should be taken into account 
in a complete calculation. As pointed out by Gross [O long ago, the [2^jv]^^''[^w]*'^^ contribution is very 
small, due to chiral cancellations. The detailed numerical study of this structure performed in ref.|12[ has 
confirmed this conclusion and shown that the crossed tern, within curly brackets, is largely dominant. It is 
due to a T^ which may be expressed as 

Tr = (a(Jo + *«oi) uu, (3) 

where the Qq^ are subthreshold coefficients [nol. In this multipolar expansion, Gqi represents the fact that a 
nucleon can be deformed when interacting with pions and is associated with the so called axial polarizability 
parameter, given by a a = 2 aj,^. 
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Figure 1: Leading contributions to the two-pion exchange potential; the black dot indicates that one of the 
nucleons acts as scalar source, that disturbs the pion cloud of the other. 



This allows the asymptotic central potential to be written as 



^a''W=[Ko+^4l) 



u u\ 



ajq(i) 



^^ 



(2) 



+ (1^2), (4) 



where crN(t) is the nucleon contribution to the scalar form factor [J. Quite generally, this form factor is 
defined in terms of the symmetry breaking Lagrangian Csb as {p'\ — Csb\p) — cr(i) u u. As discussed by 
Gasser, Sainio and Svarc [^, aN{t) is associated with the diagram of figj^a and given by 






2 m (47r)2 
where the loop integrals J{t) are written as 

[Aixf f d^Q 



JcAt) = 



JcAt) ~ 4%it) 
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J!i%it) 



(47r)2 1 



y2 



2mV-Q 



(27r)4 [(Q-q/2)2-;,2] [(g + ^/2)2_^2] [Q2 + 2TOV^.Q-g2/4] 



(5) 



(6) 
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The integral given by eq.(|g) is divergent and has to be regularized. However, this procedure affects only 
the zero-range properties of the potential and is irrelevant for the present discussion, which concentrates on 
asymptotic properties. Therefore, in the sequence, we consider only differences between functions calculated 
at generic values of t and at i = 0, denoted by a caret. 

The non-relativistic potential is obtained by going to the center of mass frame, where t = <l^ , and 
making [w u] — > 1 in eq.(0). One then has 



K^K9)--2 [a^^-q'a^^ 
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(8) 



where we now indicate explicitly the dependence of the potential on the nucleon mass. The minus sign in 
front of this expression was introduced in order to account for our convention of the relativistic T-matrix. 
It is worth noting that this contribution was already included in the general formulation of ref . |12f| . 

We define the "small" dimensionless quantities a = ^i/m, k = \q\/m, and use the standard techniques 
for loop integration in order to write 



Jc,c{^^'^^) = — da ln[l + a(l — a)K /a ] 
Jo 



2- A/l + 4a2/K2 In 






(9) 



f(i) 



^caivKg) = 2/ da{l-aY dh{l-h) 



X {[a{l-a)hK^ + {a + b-ah)a^ + {l-af{l-bf] ^ 
- [{a + b-ab)a^ + {l-af{l-hfY^\ 



Jc,c("i;9)-( — ) lN{m;q) , 
\ m J 



(10) 



^In rcf. ]18| a similar expression was used, based on the parameter ciqq = [M(agQ + 4/i'^ <^oi)l' ™tiich tends to the present one 
for very long distances. 



where the integral Iiq[m\ q) is given by 

lN{m;q) = 



da 



(a^ + a fi-?) 



tan 



tan 



■^fJ-J Jo [ \/a^ + a(l-a) k^ - (^2 + a ^2)2/4 

(1-a)- (a2+aK2)/2 \ 



I A/a2 + a(l-a) k2 _ (q,2 + a ^2)2/4 

' {o?+aK^)l2 \ 

A/a2 +a(l-a) k2 - (q;2 + a k?-Y j^j 
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tan 
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(11) 



The function crjv(i) and the asymptotic central potential can then be written simply as 

2 



K(n^;g) = - 
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[/2(a+ +ta+)] aAr(t) 



[/^ («TO - 9^ «e}i)] 



32 






lN{m;q) 



(12) 



The first result is very general and independent of specific models or approximation schemes used to 
calculate the scalar form factor. It has a rather transparent physical interpretation: the factor proportional 
to (TN(t) represents, as one will show in the sequence, the pion cloud of one of the nucleons, which is extended 
in space; the term involving the subthreshold coefficients in the square bracket represents the other nucleon, 
acting as a scalar source expanded in multipoles; the factor 2 arises because the potential is due to the crosed 
term in (TV + i?)^ and the scale (— /^^^) corresponds to the energy density of the quark condensate. 



Table 1: Values for the subthreshold coefficients Aqq, Oq]^, and the scalar form factor at the Cheng-Dashen 
point, (j{t — 2fP). Labels (emp), (s), and (A^) denote respectively empirical results, theoretical scalar form 
factor and leading delta contributions. 





emp [|10| 


emp dig 


empH 


s 


A, 


s + Ae 


ooo (m ^) 


-1.46 ±0.10 


-1.30 ±0.02 


-1.27±0.03 


0.72 


-1.89 


-1.20 


4i (m~^) 


1.14 ±0.02 


1.35 ±0.14 


1.27 ±0.03 


0.12 


0.94 


1.12 


a{2fi') (MeV) 


64 ±8 


88 ±15 


90 ±8 









Empirical values of the subthreshold coefficients can be obtained by means of dispersion relations and 
some results available in the litterature are given in table |l|. As pointed out long ago|2J, they may be well 
represented by a model including the delta and supplemented by information about the scalar form factor. 
In order to use this information, one writes a{t) = a + t a' , where a is the ttN CT-term. The parameter 
cr', the slope at the origin, may also be obtained by means of dispersion relations. For instance, a recent 
analysis |22J, based on (T(2/i^) = 60 MeV, has produced cr = 45 MeV, which corresponds to a' — 0.054^~^. 
The scalar form factor is then related to the subthreshold coefficients by Ofjpjs — a/f^ and aQi)s — o'' /f^, 
where /^r is the pion decay constant. 



As far as the delta is concerned, one uses the results of ref.l 



and obtains 
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M^ 
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(m + il/A) 



9 J Ml 






Ml 



Ml 



(13) 



+ 



2Ml{Ml 



A.2) 



AfA(3MA-TO)-^- 



2 (3Ma + 2m) 



(Ma + m) 



(14) 



where Ma is the delta mass, ^a is the pion-nucleon-delta coupling constant and we have used Z — —1/2. 
The constant ija can. be related to g by means of the large N^ result 2m,g^ = 3g/\/2 p3 and the leading 
delta contributions may be written as 



flooJA 



1 



"«i)^ = 2 Im 



fJ- 



A/a — rn 
2 1 



Ma — rn 
Adopting g — 13.40 and /^r ~ 93 MeV, one has the values displayed in table 1. 



(15) 
(16) 



Most of the dependence of the potential on the momentum transferred is embodied on the dimcnsionlcss 
function /Ar(r7T,; q), which is also influenced by the value of the nucleon mass. In the heavy baryon (HB) 
limit, defined as jti — > oo, this function can be evaluated analytically and one has |fi3 



iNim -^ oo;q) = /ffs(g) 



^) 1 ' - 







aq 



y^u^ + a{l — a) q^ 



l^'' + q^l2 



/^ m 



tan- M _ i 
2^ 2 



(17) 



The expression for the heavy baryon potential in momentum space is then 



*-"' - 1 if. 



-'oo 



2 + ' 
1 "01, 



A^^ + gV2 



tan-i ( P 



2^ 



(18) 



where we have used the Goldberger-Treiman relation g/m — gA/ fir- This result reproduces, up to irrelevant 
polynomial terms, the potential derived in ref.[|l3| using HBChPT, when one identifies the subthreshold 
coefficients and the low energy constants q as follows: a^^ = — /i^(4ci — 2c3)//^ and a^-^ = — ca//^. It 
is worth pointing out that the numerical values recommended in that work yield Oqq = —2.65 /i^^ and 
aQ2 = 1.66 /i^'^, which deviate from those quoted in table 1. This happens because the constants Ci include 
both leading terms and first order corrections p4[. The Ci then take the latter into the potential and give 
rise to fourth order corrections to otherwise third order expressions. When these corrections are eliminated 
for the sake of a consistent power counting, one recovers values compatible with table 1. 

Heavy baryon chiral perturbation theory was also used in other calculations pTl , based on time ordered 
perturbation theory. In order to allow comparison with our results, we perform Cauchy integrations in 



eqs.(||) and (0) and obtain 

lN{m\q)^ — [Jc,c~Jcs 
nu \ 



1 



2to 



{mV^ + Eq)V° + V-Q 

EQ[{mV^ + EQY-ujl][{mV^ + EqY-ujI] 

( oj+V^ + V-Q oj-V^ + V-Q 



{uol -ujI) \uj+ [(rnyo - uj+Y - E^] uj^ [{mV° ~ uj-f - E^] 



fi'^ + Q , uj- -^ yn'^ + Q 



where u}± = J ^'^ + {Q^ ± q'^/A) and Eq = y^m'^ + {Q + mV)'^. 



Going to the limit ?7i ^ oo, one has 



47r f d^Q 



iHsiq) ^ — , .„ .3 



Q" - qV^ 



Q' 



and a potential given by 



VHBig) = -l{^y{4o-eai,] 



(27r)3 



Q" - q'^/A 



Q' 



9 9 



it^' + Q 



2^2 



(19) 



(20) 



(21) 



In order to check the consistency of this result, we have integrated analytically eq.(po|) and recovered 
eq.(|l7|) exactly. 



In fig.0 we display the behaviour of ii^{m; q) for m = ttin, 10 rrijv, 100 mpf and 1000 niN, where rriN = 
938.28 MeV is the empirical nucleon mass. The ratios lN{m; q)/lHB{q) for the same values of m are shown 
in fig.^ where it is possible to see that typical deviations from the heavy baryon limit are respectively about 
30%, 10%, 2% and less than 1%. Thus, one learns that the heavy baryon limit is about 30% far from the 
real world. As we discuss in the next section, this is related to the fact that the heavy baryon expansions of 
loop integrals sometimes fail to converge m. 



The potential in coordinate space is given by 
l//(m;r) = - 



/.^M^ 



[Z'Ko + V'O] ajv(m;r). 



(22) 



where tTAr(?7i;r) is the form factor in configuration space. It is convenient to write it in terms of a dimen 
sionless function FN{m;r) as 

3^6 , „ , 2 



with 



crAr(m;r) = — ^ ( — ) FN{m;r) 
d\ 



All 
FN{m;r) = -^ 



(27r)3 



-'9-^/Ar(m;g) 



(23) 



(24) 



In the heavy baryon limit this expression can be evaluated explicitly p3| and, using x = jir, one writes 



FN{m^ 00; r) = pHsir) 



\\ ' 



3 11 

4 1 + - 



'2x 



l\ e 



X 2 x^ 



X / X 
„-2x 



(25) 
(26) 



tr 




Figure 2: Dependence of the function Iiq{m] q), eq.([lj), on the nucleon mass, for m — mjf (continuous hne), 
m = 10 rriN (dot-dashed hne), m — 100 ttin (dashed hne), m — 1000 tun (dotted hne), where tun is the 
empirical nucleon mass. 



For generic values of m, it is more convenient to perform the Fourier transform using directly eqs. 
(|l^) and one has |l| 



where 



t \' ^ m 

FN(m;r) 



1 .1 
da / db 



T^A^ Jo Jo 



X2P+2 g-A. 



2to2 ^^ (l-a)(l-6) e-''^ 



„2p 



a b 



and 



(27) 



A^ = 1/ [ail - a)b] 



2 ™2 /,,21 \2 



77^ = \a + b-ab+il-ay{l-bym^/fi'\ X 



(28) 
(29) 



The behaviour of the ratio F]\j{m; r)/ pHsif') as function of m is displayed in fig.H, where one finds the 
same qualitative pattern as in figs.0 and H. In configuration space, however, the gap between the heavy 
baryon limit and results with the empirical value of the nucleon mass varies with distance and is about 25% 
asymptotically. The difference between the sizes of the gaps in momentum and configuration spaces can 
be understood by noting that the former may contain pure polynomials, which correspond to zero-range 
functions and hence are not present in the latter. 
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Figure 3: Dependence of the ratios of the functions 1^(771; q), eg. (pi]) and Ihb{q), eg. (p^ , on the nucleon 
mass; conventions as in fig.0. 



3 Chiral Symmetry 



The behaviour of the triangle integral that determines the scalar form factor has been extensively discussed 
in a recent work by Becher and Leutwyler 0. The leading term in their chiral expansion of cr(i) can be 
obtained by neglecting terms proportional to t/A-m? in eg.(|27|) and one has 



m r^ r^ ,, (1-7^72) e-''^ 

FAm-r) ^ da db ^ ' ' ' 

T^^J- Jo Jo 



a b 



(30) 



This result is eguivalent to the dispersive representation of ref . |q| , given by 



Fi{m;r) = 



dt 



(^2 _ t/2) 



Vi 



Gi{m;t, ^) 



'^/tx I li 



with 



Ge{m; t, /i) = tan 



2m ^i - 4/i2 



(31) 



(32) 



In order to study the structure of this integral, we rewrite it as 

Fe{m; r) = - ^ ^ H dt ill-lM. G,{m; t, /.) e'^^'^^)-/'^ 



(33) 
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Figure 4: Dependence of the ratios of the functions FN{m;r), eq.(| 
mass; conventions as in fig.y. 



and Fhb{i'); eq.(P5[), on the nucleon 



For large values of x, the exponential in the integrand falls off very rapidly as a function of \/t and the 
integral tends to be dominated by a small region around t = 4/i^. On the other hand, the function G^ 
vanishes at this point and hence the net result is due to a compromise between these two effects. In the 
work of Becher and Leutwyler one learns that it is possible to write accurately 



Gi{m;t,^i) ^ 



HB 



a{t - 2^2 



2^Jt-A^i'^ 



JA 



Y^/t ^/i 



2v/i-4^2 2^ 



tan 



a/i 



Vt^^M^ 



(34) 



Js 



where a = pL/m. The first two terms in this expression correspond to the standard heavy baryon series 
whereas the last one implements the correct analytic behaviour around the point t = 4/i^. In particular, the 
HB term gives rise to eq.(p5|) and has been used in HBChPT calculations of the NN potential. 

In fig.plwe display the ratios of the functions Fhb, Fhb+Fa and Fhb+Fa+Fb, obtained by using eq.(p3) 
into eq.(^, by Fe. The first thing to be noted is that there is a gap between Fhb and Fe, which is about 
20% in the region of physical interest and can be as large as 40% asymptotically. The inclusion of F^ brings 
the approximate result much closer to the exact one, indicating that the last term in eq.(M) is important 
mainly at rather large distances. Thus, in principle, HBChPT could allow a precise numerical description 
of the NN potential, but this would require the extension of existing calculations to higher orders. This 
illustates the advantage of the relativistic framework proposed recently by Becher and Leutwyler, namely 
that it gives rise automatically to leading contributions which are already rather accurate. 
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Figure 5: Ratios of the functions Fhb (dashed line), Fhb+Fa (dotted line) and Fhb+Fa + Fb (continuous 
line) , obtained by using eq.(^) into eg. (pi]), by Fi, eqs.(|30|-p2l). 



The leading chiral contribution to the central potential is obtained by using eqs.(|3l]) and ( |32| ) into eqs.(|23|) 
and (p3). In order to assess the phenomenological implications of this result, one could evaluate phase shifts 
and compare them with experiment. However, one should bear in mind that, in contributing to phase shifts, 
the central TPEP is superimposed to both isospin dependent components and a large OPEP background. 
Its effects are completely blurred in the ^Sq wave, but may be identified in the waves ^D2 and ^G^. In this 
situation, it becomes possible for different chiral models to accomodate experimental data[|ll| |l^, |l^, Q. In 
fact, in the case of the TPEP, phase shifts do not represent the most stringent test possible. 

Therefore, in this work we compare the predictions of the asymptotic V^{Tn;r) with the Argonne phe- 
nomenological potential Q, which is rather accurate and constrained by a large data base, that also includes 
information about the ^Sq channel. In figj^ we plot the ratio of the leading chiral potential by the corre- 
sponding Argonne component, in the range lfm< r <4fm, for the sets of subthreshold coefficients considered 
in table 1. Typical discrepancies are of the order of 20%, indicating that the very simple mechanism de- 
scribed in the previous section is indeed responsible for the leading contribution to the central potential. The 
existing gap has the order of magnitude expected from chiral corrections. It is worth pointing that, in the 
case of NN interactions, these corrections would affect not only the scalar form factor, but also require the 
evaluation of other amplitudes. In particular, in the case of the delta contributions, one has also to include 
box and crossed box diagrams. A comprehensive discussion of this problem will therefore be postponed. 
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Figure 6: Ratio of the leading asymptotic relativistic chiral prediction Va{r), eqs.(|2^ 
the central medium-range Argonne potential Pq] , for the ^ 
(continuous line), ||l^ (dashed line) and Ig^l (dotted line). 



D amd (gOf^), by 
the central medium-range Argonne potential ||2^, for the values of the subthreshold coefhcients of refs. ||l^ 



4 Other Results 



We begin by considering first order corrections to the leading result for a(t), which come from two different 
sources. One of them corresponds to the terms neglected when passing from eq.(p7|) to eq.(|3n), whereas 
the other one is implemented by subthreshold coefficients, representing the remainder R discussed in the 
previous section. Thus, the scalar form factor, including both leading and first order corrections, is given by 



a{m; r) — <7N{m; r) + aR{ni; r) 



(35) 



where utv is the result of ref. |^ , given by eq. ( |29|) and ar is written in terms of the subthreshold coefficients 

3^4 ^ - 



d+ as 



with 



a^irn; r) 



Ai 



1287r2 m 



{d+, + td+,)+^-{^?-t|A)dl 



Fccir) , 



F (r)-— /'^!S_e-9-^J (a) 



(36) 



(37) 
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The main contributions to the coefhcients entering this resuh come from the cr-term and the A and one 
may also write 



o 4 2 

ar{m;r) w (7a+A{'m;r) = — -^ Fcc{r) +crA(m;r) . 



(38) 



1287r2 to/2 
The term cta is given by the diagram of fig.|l|b. Using the resuhs of ref.|2^ one finds, after neglecting terms 



CTA(m;r = -— ^1— Fa m; r , 

where Fa is obtained from the Fjv of eq. ( |30| ) by replacing rj^ with 

9^=rj^ + [(l-a)(l-6) {Ml ~ m^)/fi^] X^ . 



(39) 



(40) 



2.0 



1.5 



g 
"Id 



1.0 



0.5 



0.0 



r(fm) 



Figure 7: Ratios oi ctr eq.(pq - continuous line), a a eq.(p9|- dashed line) and (Js+a eq. (pq - dotted line) by 
the leading chiral contribution ae eqs.(p0[]3^) and ^2 



As the function 6 determines the effective mass of the two-pion system exchanged between the nucleon 
and the scalar source, this result indicates that cta has shorter range than ctat. In fig.0, we display the ratios 
of CTr, Cta and CTs+a by the leading chiral contribution CT£ and one notes that these corrections remain quite 
visible at relatively large distances. 
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The Fourier transform of the scalar form factor corresponds to a mass density, that surrounds the nucleon 
and is due, in part, to the pion cloud. As the density of the quark condensate is negative, the sign of the 
scalar density indicates that the pions, as Goldstone bosons, destroy the condensate in order to exist. In this 
picture, the energy density of the scalar cloud cannot exceed that of the condensate. In fig.^, we plot the 
function a{r)/f^fi'^ as function of the distance, and it is possible to see that this critical situation occurs close 
to 0.6 fm. On the other hand, the value 1/100, for which the pion field may be considere as weak, occurs 
around 1.5 fm, suggesting that this is the order of magnitude of the distances beyond which perturbative 
calculations may be considered as reliable. 




Figure 8: Spatial dependence of the scalar form factor, normalized by the quark condensate density: (tn 
eq.(E3 - dashed line), cta eq.pfl - dotted line) and cttv+a (continuous line). 



In order to explore this picture further, we recall that the nucleon cr-term is related to (7{r) by 

/>oo 

/ dr rV(r) . (41) 

Jo 



a = Att 



The evaluation of this integral, using <jir) — aNir) + (Ta(7') and replacing this function by /^/i^ whenever 
cr(r) > 1, yields cr = 46 MeV, quite close to the value quoted in ref.|Q. This somewhat surprising result 
is associated with the factor r^ in eq.(|34|), that suppresses the integrand at small distances, and seems to 
indicate that other mechanisms contribute little to scalar dynamics. 



Another interesting aspect of expression (|39| ) is that it allows one to study the behaviour of the scalar 
form factor in the large Nc limit, where the nucleon and the delta become both heavy and degenerate. In 
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that limit one has cta = 2(7 jv a-nd hence, 

9 fi'^ 



a{r)^a^{r)+a^{r)^^[^y 



V^ — 

X 



(42) 



This expression is identical with that obtained in the framework of the Skyrme model, where Csb{x) — 
f^^^coaF(x) and, for large distances, one hasPSJ 

\ Sttto/tt J dx X 

This result, when combined with the fact that the model accounts well for the OPEP, shows that skyrmions 
do have asymptotic properties consistent with precise perturbative calculations in the large Nc limit. On 
the other hand, in the perturbative approach, the departure from this limit is very important, suggesting 
that some of the quantitative limitations of the Skyrme model may be associated with the intrinsic iV-A 
degeneracy. 



5 Summary 

In this work we studied the central NN potential at large distances and have emphasized the special role 
played by the nucleon scalar form factor, which describes the mass density of its pion cloud. The effect of this 
cloud over the other nucleon, represented as a multipole expansion, dominates the asymptotic interaction. 
This rather simple picture accounts for about 75% of the potential, as given by the precise phcnomcnological 
fit by the Argonne group. 

We have evaluated this leading contribution, due to the exchange of two pions, and determined its 
dependence on the nucleon mass, using the relativistic formalism. Our results coincide with those of heavy 
baryon chiral perturbation theory when the nucleon mass is taken to infinity, but become about 20% smaller 
when the empirical mass is used. This shows that it is important to keep the nucleon mass finite in precise 
descriptions of the NN system. 

In a recent work, Becher and Leutwyler have shown that the scalar form factor is determined by an 
intermediate ttN amplitude, close to the threshold at t = 4/i^, which lies almost on top of a branch point. 
As a consequence, the heavy baryon expansion does not make sense for the scalar form factor. We have 
investigated the implications of this result for the NN interaction and found out that threshold effects 
can be seen at very large distances only. Therefore, in principle, HBChPT could account numerically for 
the discrepancy between the heavy baryon limit and reality found in the central potential, but this would 
require the extension of existing calculations to higher orders. This stresses the advantage of the relativistic 
framework, namely that it does not deal with truncated expressions and hence gives rise automatically to 
leading contributions which are already rather accurate. A fully covariant description of the NN interaction, 
including other degrees of freedom, will be presented elsewhere. 
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